Few-boson tunneling in a double well with spatially modulated interaction 
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We study few-boson tunneling in a one-dimensional double well with a spatially modulated interac- 
tion. The dynamics changes hom Rabi oscillations in the non-interacting case to a highly suppressed 
tunneling for intermediate coupling strengths followed by a reappearance near the fermionization 
limit. With extreme interaction inhomogeneity in the regime of strong correlations we observe tun- 
neling between the higher bands. The dynamics is explained on the basis of the few-body spectrum 
and stationary eigenstates. For higher number of particles, A?^ > 3 it is shown that the inhomogene- 
ity of the interaction can be tuned to generate tunneling resonances. Finally, a tilted double-well 
and its interplay with the interaction asymmetry is discussed. 
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I. INTRODUCTION 

Ever since the experimental realization of Bose- 
Einstein condensation, ultracold atoms have been used 
to study an enormous diversity of quantum effects with 
an unprecedented degree of control P, [1| . Advancement 
in tools Hke optical lattices Q offer us the chance to 
explore e.g. second-order tunneling Q, quantum phase 
transitions Q and non-equilibrium quantum dynamics of 
driven systems 

The double well especially serves as a prototype sys- 
tem to study interference and tunneling in great detail. 
For instance the tunneling dynamics of a Bose-Einstein 
condensate has been observed to undergo Josephson os- 
cillations 043 in which the population simply tunnels 
back and forth between the two wells. However when the 
interaction is raised beyond a critical value the atoms 
remain trapped in one wellja non-linear phenomenon 
known as self trapping 0, [13 ■ 

Of special interest are systems in lower dimen- 
sions which often display unique features. Quasi-one- 
dimensional (ID) Bose gases have been prepared exper- 
imentally by freezing the transverse degrees of freedom. 
There it is possible to tune the interaction strength be- 
tween the atoms by either using confinement induced res- 
onances or magnetic Feshbach resonances [l3 • Thus 
one can study the crossover from a weakly interacting 
to a strongly correlated regime. A particularly interest- 
ing case is the Tonks-Girardeau gas appearing in ID in 
the limit of infinitely repulsive short-ranged interaction 
which has been recently observed experimentally (l3l.[T3|. 
This gas of impenetrable boson is isomorphic with that of 
free fermions via the Bose- Fermi mapping [l5| and all the 



local properties are identical to the free fermion system. 
The gas still retains the bosonic permutation symmetry 
and so the non-local quantities differ from the fermionic 
case. 

Theoretically the quantum dynamics in the weakly in- 
teracting case has been studied using the Bose-Hubbard 
model assuming the validity of a lowest band approxima- 
tion [iGl - ligj - These studies illuminate relevant tunneling 
mechanisms and resonances. However, to capture the 
rich physics present in the stronger interaction regime 
we need to go beyond the Bose-Hubbard limit. Moreover 
numerically exact calculations of the quantum dynamics 
for few bosons through a one-dimensiona,l potential bar- 
rier (2^ or a bosonic Josephson junction (23| reveal devi- 
ations from the results obtained with mean-field calcula- 
tions as well as establish a difference between the dynam- 
ics in attractive and repulsive bosonic systems |24| . The 
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crossover from the uncorrelated to fermionization regime 
has been investigated for few bosons [2^ [13 and reveals a 
transition from Rabi-oscillations to fragmented pair tun- 
neling via a highly delayed tunneling process analogous 
to the self-trapping for condensates. The quantum dy- 
namics of an asymmetric double-wells while keeping a 
constant interaction strength has been explored in refs. 

[i3,[il,[2i,[23. 

In this investigation we go one step further and envi- 
sion a new approach to asymmetry by introducing an in- 
homogeneous, i.e., spatially varying interaction strength. 
This can be achieved experimentally by employing mag- 
netic field gradients in the vicinity of Feshbach resonances 
or by combining magnetic traps with optically induced 
Feshbach resonances (T2| . [27| . We will demonstrate how 
spatially varying interaction strengths enrich the tunnel- 
ing dynamics in the most fundamental case of a double- 
well. This has to be seen as a potential ingredient for 
more complex problems such as quantum transport in 
optical lattices. Specifically we study the crossover from 
the non-interacting to the fermionization limit for a fixed 
inhomogeneity ratio of interaction and the effect of vary- 
ing inhomogeneity ratio. An interplay of suppression and 
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resumption of tunneling is observed. For three or more 
particles, the interaction asymmetry can be tuned to gen- 
erate various many-particle tunneling resonances. Lastly 
we examine a tilted double-well and investigate the in- 
terplay between the tilt and inhomogeneity to generate 
tunneling resonances. 

The paper is organized as follows. In Section |TT] we dis- 
cuss our model and setup. In Sec. Illll we briefly describe 
our computational method. Subsequently we present and 
discuss the results for tunneling in a symmetric double 
well for two atoms first (Sec. IIV[) followed by more atom 
systems (Sec. |V]). In Sec. IVII we discuss the case of an 
asymmetric double well. 



fashion. This is computationally challenging and can be 
achieved only for few atom system. Our approach is the 
Multi-Configuration Time Dependent Hartree (MCTDH) 
method |3ll . [3^ being a wave packet dynamical tool 
known for its outstanding efficiency in high dimensional 
applications. 

The principle idea is to solve the time-dependent 
Schrodinger equation 

as an initial value problem by expanding the solution in 
terms of Hartree products <I>j = ipj ^ (g) . . . (x) (pjj^ : 



II. SETUP AND INTERACTIONS 

Our Hamiltonian (for N particles) is given by (see [sOj 
for details) 



1=1 Kj 

The double well U{x) = ^x^ + hSuj{x) is modeled as a 
harmonic potential with a central barrier shaped as a 

Gaussian 6i^(x) = ^—7= (of width w = 0.5 and height 

h = 8, where dimensionless harmonic-oscillator units are 
employed throughout). 

For ultracold atoms only the s-wave scattering is rele- 
vant and the effective interaction in ID can be written as 
a contact potential [Tlj which we sample here by a very 
narrow Gaussian. We focus on repulsive interaction only. 

The inhomogeneity of the interaction is modeled as [30| 

ff(i?)=5o[l + «tanh(f)], 

where 2R = Xi + Xj and L is the modulation length 

which we fix at L = 1. 

For R:^ L, g takes the asymptotic values 

9± = .go(l ±a). 

Thus the parameter a regulates the relative difference in 
interaction strength between the left and the right well, 

^9 = 1.9+ = 2goa, 

and the corresponding ratio is given by 



vl/(i) = ^^j(t)<i>j(t). (2) 

J 

The unknown single particle functions ipj{j = 1, 
where n refers to the total number of single particle func- 
tions used in the calculation) are in turn represented in 
a fixed primitive basis implemented on a grid. The cor- 
rect bosonic permutation symmetry is obtained by sym- 
metrization of the expansion coefficient Aj. Note that in 
the above expansion, not only are the coefficients Aj time 
dependent but also the single particle functions (pj. Us- 
ing the Dirac-Frenkel variational principle, one can derive 
the equations of motion for both Aj and <I>j. Integrating 
these differential equations of motion gives us the time 
evolution of the system via ©. This has the advantage 
that the basis ^j{t) is variationally optimal at each time 
t. Thus it can be kept relatively small, rendering the 
procedure more efficient. 

Although MCTDH is designed primarily for time de- 
pendent problems, it is also possible to compute station- 
ary states. For this purpose the relaxation method is used 
[331 . The key idea is to propagate a wave function ^0 by 
the non-unitary operator e"^"^. As r — > 00, this expo- 
nentially damps out any contribution but that stemming 
from the true ground state like £-(^"^-^0)^ ^ jj-^ practice 
one relies upon a more sophisticated scheme called the 
improved relaxation (3^ . [35| which is much more robust 
especially for excited states. Here (\E'|i/|'I') is minimized 
with respect to both the coefficients Aj and the orbitals 
(fij. The effective eigenvalue problems thus obtained are 
then solved iteratively by first solving Aj with fixed or- 
bital ipj and then optimizing ipj by propagating them in 
imaginary time over a short period. This cycle is then 
repeated. 



IV. TUNNELING DYNAMICS FOR TWO 
BOSON SYSTEM 



III. COMPUTATIONAL METHOD 

Our goal is to study the bosonic quantum dynamics 
for weak to strong interactions in a numerically exact 



We first focus on the tunneling dynamics in a sym- 
metric double-well with two bosons initially (t = 0) pre- 
pared in the left well. This is achieved by adding a tilt 
or a linear potential dx to the Hamiltonian hence mak- 
ing the left well energetically favorable. Instantaneously 
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time t 

Figure 1: (color online) Population of the right-hand well over 
time, PR{t), for different interaction strengths at q = 0.2 for 
two bosons. Inset: Long time behavior for very low interac- 
tion strength go = 0.005. Barrier height h = 8 and width 
uj = 0.5 has been used for all calculations, (all quantities are 
in dimensionless harmonic oscillator units throughout). 

the ground-state is obtained by applying the relaxation 
method (imaginary time propagation). For reasonably 
large d, this results in achieving a complete population 
imbalance between the wells. With this state as the ini- 
tial state, the tilt is instantaneously ramped down [d = 0) 
at i = to study the dynamics in a symmetric double- 
well. Our aim is to study the impact of the correlations 
between the bosons on the tunneling dynamics both with 
respect to the interaction strength as well as the spatial 
inhomogeneity. We start by fixing the inhomogeneity 
to a = 0.2 (with the left well having lower interaction 
than the right) and analyze how the dynamics varies with 
changing interaction strength go- 

A. Dynamics from the uncorrelated to the 
fermionization hmit. 

In the absence of any interaction go = 0, the bosons 
undergo Rabi oscillations between the two wells. This 
is characterized by complete tunneling of both bosons 
between the two wells with a single frequency and can be 
quantified by the time variation of the population of the 
atoms in the right well 

Pnit) = (e(.T))^(,) = / p{x;t)dx 
Jo 

where p is the one-body density. Figure [1] shows that 
Pji oscillates sinusoidally between and 1. If we in- 
troduce a very small interaction go = 0.005 (inset) the 
Rabi oscillations give way to a beat pattern due to the 
existence of two very close frequencies. Increasing the 
interaction strength further (go = 0.2), we observe a sup- 
pression of tunneling with the maximum population in 
the right well p^™"^ « 0.2. This is a manifestation of 
the inhomogeneous interaction which drives the tunnel- 
ing off-resonance and should be carefully distinguished 
from the delayed pair-tunneling and self-trapping for the 
same go-value in the case of homogeneous interactions 



(see below). The dynamics consists of a slow tunneling 
envelope with suppressed amplitude, which is modulated 
by a faster oscillation. For higher values of interaction 
strength (go = 4.7), the tunneling is completely sup- 
pressed. What remains is a fast oscillation with a tiny 
amplitude. 

However, contrary to the naive intuition a reappear- 
ance of tunneling occurs for larger values of the coupling 
strength. We observe a partial restoration of tunneling 
with p^™"^ = 0.7 for the value go = 150, which is close 
to the so called fermionization limit. The dynamics is 
characterized by two frequencies - one very close to the 
Rabi frequency modulated by a faster oscillation. Ideally 
at the fermionization limit go — > oo, the system of hard- 
core bosons maps to a system of free fermions [l^ and all 
the local properties are identical. Hence in this limit we 
would have complete two-mode single particle tunneling 
analogous to tunneling of two free fermions. 

Before we move on to analyze in detail the above ob- 
servations, let us comment briefly on the differences be- 
tween the behavior observed in our setup and the case of 
a symmetric double-well with homogeneous interaction. 
Clearly both for the non- and infinitely interacting lim- 
its the inhomogeneity doesn't play a role. For homoge- 
neous interactions and a symmetric trap tunneling is al- 
ways resonant and complete. However, different strength 
of interaction yield different dynamics like a transition 
from pair-tunneling for low interaction strength to a 
self-trapping mechanism for larger interaction strength 
which is characterized by extremely long tunneling times 
[l6l . [20I . [2ll . [25l . [26| . In our case though we observe an 
actual suppression of the tunneling amplitude and not so 
much a delayed process. In case of an asymmetric-well 
with homogeneous interaction, the effects in the low in- 
teraction regime are equivalent to our set-up: The tilt 
has the same effect as an interaction asymmetry, namely 
it destroys resonant behavior thereby leading to a sup- 
pression of tunneling (TtI. [18|. Nevertheless, our case is 
fundamentally different and this is evident in the strong 
interaction regime. Specifically the reemergence of tun- 
neling we observe does not occur in the tilted double- well 
system. 



B. Analysis 

The understanding of the above-described dynamics 
lies in the variation of the few body spectrum as go is 
changed from zero to the fermionization limit (Figj^l^a)). 
Considering the wavefunction 'i!{t) = X^m 
with energy E„i corresponding to the stationary state 
the population imbalance S{t) = (0(a;) — 0(— a;)),p(t) 
can be computed to be 

cos(a;m„ 

m<in m 
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Figure 2: (color online): (a) Two particle energy spectrum as 
a function of the interaction strength go for a = 0.2. Inset: 
Lowest energy levels for low interaction strength. Bottom. 
Few body energy spectrum with go for (b) a = and (c) a = 
0.2. 



where W„ 



(1'™|e(a:;)|1'„)c„ 



and ( 



The energy spectrum of both the non-interacting and 
the fermionization limit can be understood from the sin- 
gle particle energy spectrum of the double well, which is 
in the form of bands each pertaining to a pair of sym- 
metric and antisymmetric orbitals. 

In the uncorrelated limit (go 0) the low-lying en- 
ergies of the spectrum are obtained by distributing the 
atoms over the symmetric and antisymmetric single par- 
ticle orbitals in the first band. This leads to TV + 1 energy 
levels, N being the number of bosons. = i?o + mA° 
with m = 0, N where A*^ = ei — eo is the energy dif- 
ference between the two single particle orbitals in the 
first band. Thus for go = the levels are equidis- 
tant (Figl2l[a) inset) and we see Rabi oscillation with 
frequency woi = '-^12 = the interaction is in- 

creased (go = 0.005), this equidistance is slightly broken 
(woi — W12) and we get a superposition of two very close 
frequencies. This results in the formation of the beat 
pattern seen in the dynamics for go = 0.005. 

To understand the dynamics in the low interaction 
regime, it is instructive to map our system to a two-site 
Bose-Hubbard Hamiltonian (2a, u3 



H 



^ 2 

= L.R 



(h, - 1) (4) 



where J is the tunneling coupling, Ul.b. is the on-site 



'3^3- 



energy of the left / right well and iij = c 

Before proceeding, we note here that there is no direct 
connection between the time-dependent SPF used within 
our numerical M.C.T.D.H. calculations and the param- 
eters of the Bose-Hubbard (B-H) Hamiltonian. In the 



standard B-H model (which is valid in the weak interac- 
tion regime) the parameters J, Ul, Ur are time indepen- 
dent constants while the shape of the orbitals such as the 
localized Wannier functions retain the shape throughout 
the course of the dynamics. Moreover, even for low en- 
ergies the two most occupied modes for propagation do 
not necessarily coincide with the two modes of the B-H 
model. In this weak interaction regime, the B-H model is 
just a good approximation to our exact calculation and 
thus we have used it as solely an explanatory tool to an- 
alyze the results. 

Using the B-H Hamiltonian for Ul^Ur ^ J, the 
highest two eigenvalues are approximately Ur and Ul- 
Whereas in the homogeneous case a = these two lev- 
els are close to degenerate Ul ~ Ur (Fig[2I^b)), here we 
have a breaking of the parity symmetry since Ur > Ul 
(FiglHc)). This is understandable since two particles lo- 
calized in the left well have lower energy than two parti- 
cle in the right well leading to the energy level separation 
seen in Figl^Jc). In terms of the number-state represen- 
tation in the localized basis \Nl^^\ Nr^^^) the degenerate 
eigenstates for the homogeneous case read 

0i,2«;)5(|O,2)±|2,O)) 

and consequently the dynamics consists of shuffling the 
probability between the two states corresponding to a 
complete two particle tunneling. 

In the case of sufficiently strong inhomogeneous inter- 
action, the removal of the degeneracy of the energy lev- 
els leads to a decoupling of the eigenstates into localized 
number-states 

(/.i«|2,0) ,</)2«|0,2) 

This implies that the initial state ?/'(t = 0) = |2,0) is 
very close to the 1st excited state and thus is effec- 
tively a stationary state of the system. This results in 
the suppression of tunneling for corresponding values of 
.90 

In the fermionization limit (go 00) the system pos- 
sesses the same local properties as a system of non- 
interacting fermions due to the Bose- Fermi mapping [l5| . 
Thus in an ideal case the inhomogeneity doesn't manifest 
(g-j- 00) and the tunneling dynamics is identical to a 
system of free fermions. As an idealization if we consider 
the initial state as two non-interacting fermions in the 
left well, then they would occupy the lowest two orbitals 
localized in the left well. In terms of the single particle 



eigenstates of the double well \n[fi^) where ria^J denotes 
the occupation number of the symmetric {ap ~ 0) or an- 
tisymmetric (ttfj = 1) orbital in band /3, the tunneling 
frequencies w„„' = En — E'^ are given by [2^ 



,(/3) 



(5) 



=o,±i 



where denotes the energy splitting of the band (3 , ni^ 
represents the occupation of the anti-symmetric orbital 
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time t 

Figure 3: (color online): Population of the right-well over 
time, PR{t), at go — 0.2 for different a values. Inset: Varia- 
tion of maximum population of the right well Pj^"^"-^ with a 
for go = 0.2. 



of the band /3. Thus for two particles the contributing 
frequencies are the lowest band Rabi frequency A*' and 
the tunnel splitting of the first excited band . The tun- 
neling dynamics can be pictured roughly as two fermions 
tunneling independently in the first two bands. 

In our system however the finiteness of the 170 value 
leads to deviations from the ideal fermionic dynamics. 
The inhomogeneity of the interaction still manifests lead- 
ing to a difference w.r.t the localized two-particle energy 
level in each well and the tunneling remains incomplete. 



C. Dynamics with varying inhomogeneity 

Having analyzed how the dynamics varies with chang- 
ing interaction strength at a fixed interaction asymmetry, 
it is worthwhile to study the dependence of the tunneling 
dynamics on the strength of the inhomogeneity. For this 
we study the effect of different a values on the tunneling 
dynamics for a fixed go ~ 0.2. 

In Fig. [3] we observe that for a = 0, we have complete 
tunneling with a two mode dynamics i.e. fast oscillations 
(luqi) which modulate slower tunneling oscillations (W12). 
When a is increased to a value of 0.04, the tunneling max- 
imum is reduced to roughly 0.7 while still retaining the 
two-mode character. As a is further increased to 0.2 the 
tunneling is suppressed as described in Sec. IIVBI The 
characteristic display of fast and slow oscillations aris- 
ing due to the time-scale difference of the contributing 
frequencies is not prominent here and for higher interac- 
tion asymmetry (a = 0.5) we have effectively single mode 
tunneling with frequency wqi. 

The variation of the maximum population with 
the inhomogeneity a (Fig[3] inset) shows a sharp drop 
with increasing a before effectively reaching a constant 
value ~ 0.12 for a > 0.3. The reader should note that 
p^max jjQgg jjQ^ gQ 2ero in the asymptotic limit a — !> 1 
or ^fi- ^ 00. This is due to the fact that with a fi- 
nite value of go and a finite barrier height the tunneling 




Figure 4: (color online): (a) Population variation with time 
Pj{(t) at go = 25 and a = 1 for Pr{0) = 1, i.e initially 
populating the right-well, (b) Energy spectrum for a = 1 



coupling (J) is not negligible compared to Ur. As a 
consequence there remains a finite probability of bosonic 
tunneling between the two wells. 

D. Strong interaction inhomogeneity 

An extremely strong inhomogeneity at a high interac- 
tion value leads to an interesting higher band tunneling 
dynamics. We can realize such a system by having a = 1 
at go = 25. This set up effectively makes the bosons 
fermionized in the right-well and almost non-interacting 
in the left. Preparing the initial set-up with both bosons 
in the left well leads to the suppression of tunneling. 
However if we prepare the initial state with two boson 
in the right well, then we observe substantial tunneling. 
In FiglH (a) we see that the Pr oscillates between 1 and 
0.5 indicating a single boson tunneling with a single dom- 
inant frequency. 

In order to understand the phenomenon we look at the 
energy spectrum at a = 1 (Fig|4](b)). While the ground- 
state remains unaffected, what we see is that close to 
the fermionization regime (go — 25), the first excited 
state decouples from the higher three states which come 
closer. The main contribution to the first excited state 
is the state |2, 0) and its separation from the other states 
could be understood from the fact that two boson in the 
left-well is almost non-interacting and thus energetically 
far off resonant from two effectively fermionized boson in 
the right-well |0, 2). The consequences of this fact are the 
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Figure 5: (color online) (a) Population of the right-hand well 
over time, PR(t), for three bosons for different interaction 
strengths at a = 0.2. (b) Three boson energy spectrum at 
a = 0.2. 

following: (i) The initial configuration of |2, 0) becomes a 
stationary-state resulting in a highly suppressed tunnel- 
ing, and (ii) the state |0, 2) of the lowest band becomes 
energetically resonant and couples to the states 
and jl*^,!^) in the higher bands (where the superscript 
refer to the ground (0) or excited (1) orbital of the corre- 
sponding well). The latter leads to a tunneling dynam- 
ics in the higher band states predominantly between the 
2nd and the 4th excited eigenstates (see Fig. 111(b)) which 
have greater overlap with the initial state |0, 2). These or- 
bitals have mostly contributions from the states |0, 2) and 
1°) while the other orbital has minimal overlap with 
the initial state. As a result we get a single-particle tun- 
neling with one dominant frequency given by the splitting 
of the energy between these two levels. In other words, 
we effectively have a single boson tunneling between the 
wells in the excited band. Note that this highly corre- 
lated single-particle tunneling scenario is attributed to 
the high inhomogeneity in the strong interaction regime 
since the combination of these two factors are responsible 
for turning the pair-tunneling scenario off-resonance. 



V. MULTI-PARTICLE DYNAMICS 

Having analyzed the tunneling dynamics of two atoms 
let us now focus on the case of three or more atoms to 
see the general atom number dependence of tunneling in 



the presence of spatially modulated interactions. 



A. General behavior and mechanisms 

Like in the two boson case we start with the initial state 
of iV = 3 bosons prepared in the left well. As shown in 
Fig. [5]Ja), the main effects are similar to the two-atom 
case. The dynamics is again governed by frequencies de- 
termined by the energy difference of the low lying spec- 
trum. For very small interaction, the nearly equal energy 
difference gives rise to the beat pattern similar to that 
of two particles. As we increase the interaction strength, 
we observe suppression of tunneling for go = 0.2 followed 
by a partial restoration at = 4.7 and a higher am- 
plitude reemergence close to the fermionization limit at 
(7o = 150. The general mechanism for the suppression is 
the same as for the two particle case. Now, however, in 
the symmetric case a = 0, the contributing nearly degen- 
erate eigenstates are of the form | A^, 0) ± |0, A^). Conse- 
quently we have a complete TV particle tunneling with a 
frequency given by [Hi w ~ 2NU/{N - 1)! x {2A°/U)^ 
where U = Ul^Ur denotes the on-site interaction en- 
ergy. The tunnel period thus grows exponentially with 
A". When the inhomogeneous interaction is introduced, 
the states decouple to the localized number-states | A^, 0) 
and |0, A^) and thus the initial state becomes a station- 
ary one leading to the suppression of tunneling. The 
important thing to note is that with increasing A^, the 
suppression of tunneling occurs for much smaller values 
of go. For instance at go = 0.2 for A^ = 3 we have al- 
most complete suppression in contrast with N = 2 where 
we still observed significant tunneling (see FigH]) for this 
value of go- This could be understood from the fact that 
the contribution of the on-site energy on the cat-state 
goes as ~ Ul,rN{N — l)/2, while that of the tunneling 
term is A^ independent. This fact is responsible for a 
significant decoupling of these states at a lower go value 
leading to faster suppression of tunneling as A^ increases. 

Also unlike that of the two boson case, the spec- 
trum for the three boson case contains crossings between 
the higher- lying states (see FiglSJb)) and in the vicin- 
ity of these crossings there is a partial reappearance of 
tunneling. This can be seen for instance at go = 4.7 
where we observe a restoration in the three-particle case 
whereas for two particles we still observed a significant 
suppression (see Fig. [1}. In this regime the higher 
bands contribute more significantly leading to the con- 
voluted dynamics observed. These higher band contri- 
butions leads to further recovery with increasing interac- 
tion strength towards the fermionization regime although 
even for go ~ 150 we do not get the exact fermionic dy- 
namics which is characterized by the tunneling of three 
independent fermions. 
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Figure 6: (color online): Population of the right-well over 
time, Pii{t), at go = 0.2 for different a values for (a) 3- 
particles and (b) 4-particles. Inset: Variation of maximum 
population of the right well p^*""^ with a for go = 0.2. 



B. Generating tunneling resonances by interaction 
inhomogeneity 

A very interesting phenomenon for the > 3 particle 
case is that by tuning the asymmetry a we get a control- 
lable reemergence of tunneling. To observe this, we study 
how the tunneling dynamics changes with different values 
of a for go ~ 0.2 (FiglHl). The value of 50 is chosen such 
that the inhomogeneity effect manifest but is still in the 
two-mode regime. For three atoms we observe (FigEI^a)) 
that a complete tunneling for a = gives way to sup- 
pressed tunneling with increasing a value. However at 
a = 0.5 we observe a reappearance which is in form of a 
tunneling resonance peaked at a = 0.5 with Pjj™"^ 0.6 
corresponding to effective two boson tunneling. In the 
case of = 4 we see two resonances (figlHl^b) mset) - 
the larger one centered on a = 0.3333 with an amplitude 
0.75 and the smaller one at a = 0.6667 with an amplitude 
0.5 resulting in the reappearance of tunneling shown in 
FiglUb). 

In order to understand this we have to study the 
spectra and the underlying eigenstates for different a 
(Figl?]). In the case of A^ = 3 for no asymmetry a = 0, 
the highest two levels form a doublet (FiglTl^a)) and 
the corresponding eigenstates are degenerate of the form 
-^(|3,0} ± |0,3)). As a is increased the parity symme- 
try is broken and the doublets separate and likewise the 
eigenstates decouple (FigjTlJb)). The energy eigenvalues 
(in the limit of very high go) are given by Ul, Ur, 3Ul 
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Figure 7: (color online) :Three particle energy levels for < 
go < 0.3 for (a) a = 0, (b) q = 0.04 and (c) a = 0.5. 



and 3Uii with the corresponding eigenstates |2, 1), |1, 2), 
|3, 0) and |0, 3). However, when Ur « 3Ul {a = 0.5) the 
1st and the 2nd excited eigenstates become near degen- 
erate and form a doublet of the form -^(|1,2) ± |3, 0)) 

(FiglT^c)). Thus the initial state |3, 0) is no longer a sta- 
tionary state of the system. As a consequence we get a 
restoration of tunneling and the dynamics basically in- 
volves shuffling atoms between these two number-states. 
In other words we have tunneling of two particles be- 
tween the two wells while one particle remains in the left 
well. This resonant two particle tunneling is what we ob- 
serve for the a = 0.5 case. As a is increased further this 
degeneracy is once again broken and the states decou- 
ple leading back to the suppressed tunneling dynamics. 
This is reminiscent of what happens in the asymmetric 
double- well for homogeneous interactions (T7| . 

In similar consideration, for the 4-particle case the en- 
ergy eigenvalues are 3Ul, QUl, {Ul + Ur), 3Ur and 
6Ur. Now if Ur 2Ul {a = 0.3333) then we have 
two degeneracies viz 3Ur 6Ul and {Ul + Ur) 3Ul 
corresponding to the eigenstates -^(14,0) ± |1,3)) and 

^(|3,1) ± |2,2)). Since the initial state is |4,0) only 
the first degeneracy contributes. Thus the dynamics in 
this case consists of tunneling of three bosons between 
the wells while one boson remains in the left well. This 
results in the tunneling amplitude of 0.75. The second 
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Figure 8: (color online) Temporal evolution of (a) pair- 
Probability and (b) three particle probability at q = and 
a = 0.5 for = 3 and go = 0.2 . 



Figure 9: Variation of (a) tunneling maximum p^""^^ with 
tilt d (b) maximum single particle probability pi with tilt d 
for AT = 2, go = 0.2 and a = 0.2. 



tunneling peak occurs for Uu — > 5Ul {a — 0.6667) which 
leads to {Ul + Ur) — !> 6Ul- The corresponding degen- 
erate eigenstates are --i=(|4,0) ± |2,2)) and we observe 
tunneling of two bosons on top of others remaining in 
the left well and thus the tunneling peak of 0.5. The 
above analysis can be extended generically for N parti- 
cles where we would have N — 2 resonances corresponding 
to the degeneracies between the eigenstates. 

C. Correlations 

In order to study the exact nature of tunneling dy- 
namics, we need to investigate the correlations between 
the particles. For this we study the temporal evolution 
of the pair-probability or the probability of finding two 
particles in the same well defined by 



the dynamics by the eigenstate analysis in the preceding 
section as tunneling between |3,0) and |0,3) states. 

Similarly at resonance (a = 0.5) we find that os- 
cillates from 0.1 and 1 implying that the system os- 
cillates between a three-particle state to a non-three- 
particle state, namely the pair-state |1,2) which can be 
inferred from the variation of p2 (Fig[8jb)). As a result 
we have pair tunneling on top of a particle remaining in 
the left- well. (Ideally in the case of B-H model, p2 should 
be oscillating between 1 and 0.33 while p^ between 1 and 
0. However in our case the realistic potential and param- 
eter regimes as well as some higher band contributions 
leads to the some deviations from this behavior). 



P2it) = (e(xi)e(x2) + ei-xi)e{-x2))t (6) 

and the three-particle-probability or the probability of 
finding all three particles in the same well defined by 

P3{t) = (e(xi)e(x2)e(x3) + e(-xi)e(-x2)e(-x3))t 

(7) 

In the case of iV = 3, for homogeneous interaction a = 
at go = 0.2 both p2 and ps oscillate close to unity 
(FiglHI). This implies that all the three particles can be 
found in the same well or in other words they tunnel 
together between the wells. This confirms the analysis of 



VI. ASYMMETRIC DOUBLE- WELL 

Thus far we have investigated the dynamics in symmet- 
ric double- well with inhomogeneously interacting bosons. 
An interesting extension is to study the dynamics in an 
asymmetric double-well. This gives us the chance to ex- 
amine the interplay between the interaction inhomogene- 
ity and the tilt. A special interesting consideration would 
be to see if the tilt could be tuned to offset the inho- 
mogeneity in the interaction and mimic the dynamics of 
symmetric interaction case or further if it can generate 
some new tunneling resonances. 
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A. Generating tunneling resonances by a tilt. 

In symmetric wells with homogeneous interaction, the 
localized N particle state |A^, 0) has the same energy as 
that of the state |0, A^) resulting in a complete iV-particle 
tunneling between the wells. With the introduction of 
the inhomogeneity w.r.t the interaction, this resonance 
is broken and the energy of TV particles in the right well 
is higher than that in the left well resulting in the sup- 
pression of tunneling as seen before. Now if we incorpo- 
rate a tilt in the double well such that the left well is 
lifted and right well is pushed down energetically in ex- 
actly the right amount to make the localized N particle 
energy levels resonant then we should expect a reemer- 
gence of tunneling. 

To observe this we prepare the initial state with both 
particles in the right well tp{Q) = |0,2) and study the 
variation of the tunneling maximum with a tilt d 

(FigMa)) incorporated into the Hamiltonian as a lin- 
ear term —dx. We restrict ourselves to the a = 0.2 
and go = 0.2 case. We observe a sharp resonance at 
d « 0.0065 corresponding to the tilt which exactly bal- 
ances the localized pair-state energy difference due to in- 
homogeneous interaction. The result is pair-tunneling 
between the two wells as we would have it in a completely 
symmetric set-up. 

With higher tilt the tunneling maximum falls off very 
sharply as the pair-state becomes off-resonant again and 
we get a suppression of tunneling. The next maximum 
occurs when the tilt is large enough to make the local- 
ized pair state |0, 2) resonant with the state |1, 1). This 
results in a broad tunneling maximum at d « 0.045 cor- 
responding to single-particle tunneling. 

To confirm our analysis of the tunneling mechanism 
we look at the variation of maximum single particle 
probability pi with tilt (Fig[9fb)), defined as pi = 
maxi(l — j52(i)) which gives the probability of having 
only one particle in a well. We observe a negligible value 
at the first resonance d « 0.0065 confirming that the dy- 
namics is pair-tunneling while a very broad maximum 
peaked at the second resonance d « 0.045 corresponds 
to the maximum probability of finding a single particle 
which in our case is the |1, 1) state and the dynamics is 
a single particle tunneling between the |0,2) and 
states. 



B. Spectral Analysis 

To understand the effect of the tilt on the tunneling 
dynamics we study the energy spectra E with varying 
tih d at fixed go = 0.2 and a = 0.2 (FigdHl). At d = 
the eigenstates are basically number-states in the local- 
ized basis. With increasing d, the highest two levels |0, 2) 
and 1 2, 0) move closer and form a sharp avoided crossing 
at d ~ 0.0065 corresponding to the first tunneling res- 
onance. At this point the tilt exactly balances the in- 
teraction inhomogeneity and the eigenstate is in form of 
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d 

Figure 10: (color online) Two particle energy spectrum with 
tilt d for Q = 0.2 and go = 0.2. 

the cat-state |2, 0) ± |0, 2). This state is very sensitive to 
the tilt and a minute perturbation decouples them into 
the localized number-state resulting in a very sharp tun- 
neling resonance. The ground-state, which is the |1,1) 
state is insensitive to the tilt since this lowering of one 
particle and raising another particle keeps the state ener- 
getically unaffected within the linear regime. This state 
forms a broad (anti)crossing with the lower excitedstate 
at d « 0.045 forming the broad single-particle tunnel- 
ing resonance seen in the dynamics. This behavior seen 
in the two-particle case can be expected in general for 
N particles giving N resonances corresponding to the 
avoided crossings encountered. In particular with in- 
creasing tilt, the successive resonances corresponds to a 
mechanism where one less particle tunnels compared to 
that of the previous one while the width of the resonances 
becomes progressively broader. 



VII. CONCLUSION AND OUTLOOK 

We have investigated the double-well tunneling dy- 
namics with inhomogeneous interaction. More specifi- 
cally we modeled the system such that we have two dif- 
ferent interaction strengths in the two wells. What we 
observe is that this inhomogeneity leads to a suppres- 
sion of tunneling. The reason for this suppression can 
be attributed to the breaking up of the doublet struc- 
ture in the energy spectrum leading to a decoupling of 
the eigenstates into the localized number-state. Increas- 
ing the interaction to the fermionization limit leads to a 
reappearance of the tunneling. The dynamics is governed 
by the band splitting of the first two bands although the 
finiteness of the interaction strength and the presence of 
the interaction inhomogeneity leads to deviation from the 
ideal fermionic behavior. For a very pronounced interac- 
tion inhomogeneity for strong interactions, we observe 
single particle tunneling between the localized excited 
bands of the double- well. 

These basic considerations can be used to understand 
the many particle system. There we observed a more 
severe suppression of tunneling for even lower go values. 
Most importantly for > 3 atoms, one can generate 
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tunneling resonances by tuning the interaction asymme- 
try. These resonances occur as a result of the forma- 
tion of degeneracies between different eigenstates. For 
three particles, the exact tunneling mechanism was in- 
vestigated using the evolution of the pair-probability and 
the three-particle probability. These studies show that 
we get correlated pair and triplet tunneling with a com- 
plete absence of single particle tunneling. 

Finally we explored the dynamics in a asymmetric 
double-well and this gives us an understanding of the 
interplay between the interaction inhomogeneity and the 
tilt. We observe that the tilt can be tuned to offset the 
interaction inhomogeneity leading to a tunneling reso- 
nance. These dynamics have been explained through the 
spectral analysis in terms of avoided crossings between 
the levels. Note that an interesting prospective would 
be to try to describe the presently found effects in the 
context of a generalized Bose-Hubbard model, where the 
on-site energies and the coupling constants would be site- 
and occupation number dependent [s^-fs^. 

Understanding the few-body mechanisms of tunneling 



with spatially modulated interactions can be used to de- 
sign schemes for selective transport of particles between 
different wells and/or reservoir systems (39l . |40| . Fur- 
ther our study could serve as a starting point for the 
investigation of the quantum dynamics in the presence 
of time-dependent interaction modulations and even be 
extended to multi-well systems (4l| . 
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